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We assess heat and mass transfer limitations in in situ studies of model catalysts with a first- 
principles based multiscale modeling approach that integrates a detailed description of the surface 
reaction chemistry and the macro-scale flow structures. Using the CO oxidation at RuO2(110) as a 
prototypical example we demonstrate that factors like a suppressed heat conduction at the backside 
of the thin single-crystal, and the build-up of a product boundary layer above the flat-faced surface 
play a signiflcant role. 



I. INTRODUCTION 

During practical operation heterogeneous catalysts are 
exposed to rather harsh environments with reactant par- 
tial pressures of the order of atmospheres and quite el- 
evated temperatures. Quite in contrast to these con- 
ditions, much of our atomic-scale understanding of the 
catalytic function derives to date from controlled experi- 
ments of single-crystals in ultra- high vacuum (UHV) . In 
order to bridge particularly the resulting pressure gap 
much progress has recently been achieved in pushing in- 
situ techniques that are capable of delivering equally re- 
solved and quantitative information for such model cat- 
alysts in technologically relevant gas-phasesi. The fo- 
cus in such studies is on identifying possible differences 
in the surface chemistry at corresponding near-ambient 
pressures and elevated temperatures. However, one also 
has to recognize that at then achievable, much higher 
conversion rates heat and mass transfer effects in the gas- 
surface system become increasingly important. In con- 
trast to low-conversion operation in UHV this concerns 
for instance the gas-phase transport of formed products 
away from the active surface, and how efficiently the large 
amount of heat generated by the exothermic surface re- 
actions can dissipate in the system. 

Differences with respect to these factors in present ex- 
perimental in-situ setups may possibly account for some 
of the existing controversies in the field. Similarly, cor- 
responding effects are important for comparisons with 
quantitative data from emerging first-principles based 
microkinetic modeling approaches?. Among these, lat- 
est first- principles kinetic Monte Carlo (Ip-kMC) mod- 
eling even provides a microscopically correct account of 
the spatial arrangement and interactions of the adsorbed 
chemicals^. The thereby achieved predictive-quality de- 
scription of the detailed surface reaction chemistry must 
then be complemented by an appropriate treatment of 
the macroscopic heat and mass transfer in the gas-surface 
system. Here we aim to qualify the latter effects precisely 
for the near-ambient, high-activity conditions targeted by 
the novel in-situ techniques. We therefore present an ap- 
proach which efficiently integrates accurate Ip-kMC sim- 
ulations for the surface kinetics into computational fluid 
dynamics (CFD) simulations. Using the CO oxidation 



at RuO2(110) as a prototypical example, we illustrate 
that the peculiarities of the thin single-crystal reactor 
geometry can readily lead to heat dissipation and mass 
transport limitations that severely affect the observable 
catalytic function. If these limitations are not appropri- 
ately accounted for in both experiment and theory, wrong 
mechanistic conclusions about the atomic-scale surface 
kinetics at technologically relevant gas-phase conditions 
may easily arise, thereby blocking the envisioned route 
towards a rational design of future improved catalysts. 



II. METHODS 

At the continuum level a reliable description of the 
macro-scale flow structures is achieved by the transient 
Navier-Stokes equations together with energy and species 
governing equations*. Aiming for a general assessment of 
heat and mass transfer limitations we confine the atten- 
tion to the center of a flat-faced single-crystal model cat- 
alyst of thickness d and consider the axisymmetric stag- 
nation flow geometry shown in Fig. [TJ In this geometrj*^, 
the equations are closed by appropriate boundary condi- 
tions at the inlet at a macroscopic distance z = L and at 
the solid surface at z = 0. For the inlet this is obviously 
the control of the temperature, density and composition 



of the gas mixture, here T'"\ p^]^ , p^Q and Pc02' '"^ 
well as its axial velocity m'"'. With thus defined inlet 
conditions, we here specifically target stationary opera- 
tion, i.e. the steady-state catalytic conversion rates, so- 
called turnover frequencies (TOFs, in units of product 
molecules per surface active site and time). In this case, 
the partial mass fiuxes of reactants and products at the 
gas-solid interface cancel to yield a total mass flux of zero. 
For near-ambient pressures we can furthermore assume 
the temperature to be continuous across the interface. 
The surface boundary conditions that then remain to be 
determined are the species conversion rate and connected 
heat release through the on-going surface chemical reac- 
tions, as well as the heat flux into the solid. 

For the prior two quantities we rely on the accurate 
first-principles description provided by Ip-kMC simula- 
tions, which evaluate the statistical interplay between the 
elementary processes involved in the catalytic cycled. In 
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FIG. 1: Comparison of intrinsic turnover frequencies (TOFs, 
black dashed line) with observable TOFs (red solid line) when 
accounting for heat and mass transfer effects. In the consid- 
ered stagnation flow geometry the inlet is placed at distance 
L = 1 cm above the RuO2(110) surface as schematically ex- 
plained in the inset. Shown is steady-state data for T'"' — 
500K, = 0.3 atm, it'"' — 1 cm/s and varying Pqq at the 
inlet. Due to the suppressed heat flux at the back of the d — 
1mm thin single-crystal the system is able to sustain a high- 
activity branch for more CO-rich conditions than the nominal 
most active "state", where the surface temperature, shown 
in the lower panel, is significantly increased. 



the established model of CO oxidation at RuO2(110) this 
is specifically the set of 26 elementary processes defined 
by all non-correlated site and element specific adsorp- 
tion, desorption, diffusion and reaction events that can 
occur on a lattice spanned by two different active sites 
offered by the surface^i^. For a given gas-phase impinge- 
ment as described by a local temperature and reac- 
tant partial pressures Pq^ and p^Q directly at the surface, 
these simulations provide the average rate of reactant and 
product adsorption, desorption and conversion per sur- 
face area and time, i.e. the partial mass fluxes at the 
surface required for the CFD modeling. Multiplying the 
conversion rate with the specific enthalpy change con- 
nected with one CO + O ^ CO2 reaction event provides 
furthermore the released heat rate that enters into the 
heat balancing equation. While it is only these average 
fluxes and heat rate that matter for the macroscopically 
described flow field, it is important to note that they are 



still properly derived from microscopic simulations that 
fully account for the site heterogeneity and distributions 
at the surface. This is thus distinctly different to mean- 
field based phenomenological descriptions that are com- 
monly integrated in the CFD modelling of macro-scale 
flow structures^, and which have for the CO oxidation at 
RuO2(I10) been shown to fail qualitativelyi. 

The also required heat flux into the solid is determined 
by the heat transport inside the crystal of width d (de- 
scribed through the bulk heat conductivity) and the de- 
gree of heat dissipation that is possible at the back of 
the crystal, e.g. through radiative loss or contact with 
the sample holder. With real experimental setups ly- 
ing anywhere in between we analyse the relevance of this 
factor for thin single-crystals by considering two opposite 
extremes: Fixed temperature at the sample backside to 
mimic a highly efficient heat coupling of the crystal to 
the system, and zero heat flux at the sample backside to 
represent a well insulated sample. 



III. RESULTS AND DISCUSSION 

In the adiabatic situation with a suppressed heat flux 
through the back of the sample the main dissipation 
channel left for the heat released by the exothermic sur- 
face reactions is into the surrounding gas-phase. Above a 
critical value of the TOF (and therewith generated heat 
rate) this channel may no longer be efficient enough to 
maintain the nominal surface temperature. If increas- 
ing the surface temperature then furthermore enables 
higher conversion rates, the system can run away into a 
new highly active steady-state that is characterized by 
sizable temperature gradients across the system. Con- 
sidering representative parameters for the crystal width 
d = 1 mm, inlet distance L = 1 cm and axial inlet flow 
velocity u'"' — 1 cm/s in our simulations for the CO ox- 
idation at RuO2(110) we find this critical TOF to be 
of the order of 10 site~^ sec^^ at near-ambient partial 
pressures. Variations of d, L and u'"' over one order of 
magnitude change this value by a similar amount, but 
leave the qualitative physics discussed in the following 
untouched. 

In the most active "state" intrinsic steady-state con- 
versions above this critical TOF number are reached for 
temperatures above about 500 K, and Fig. [T] illustrates 
the sizable effect that the surface heat-up already has 
on the observable activity at this threshold. Compared 
are the intrinsic steady-state TOFs as resulting from the 
Ip-kMC simulations, i.e. assuming that the tempera- 
ture and partial pressures at the surface are identical 
to those at the inlet, and the really observable steady- 
state TOFs when explicitly accounting for heat and mass 
transfer effects through the coupled lp-kMC-|-CFD ap- 
proach. Shown is a set of gas-phase conditions for fixed 
Pol ~ ^■'^ ^^'^ that comprises with increasing p™Q the 
three characteristic "states" of the surface: 0-poisoned 
at the lowest p'^Q shown, CO-poisoned at the highest 



Pqq shown, and in between the most active "state" with 
both reactants coexisting in appreciable amounts at the 
surface^. For a range of more CO-rich conditions than 
this intrinsically most active "state" of the surface, the 
Ip-kMC+CFD results in Fig. [T] reveal that the system 
is in fact able to sustain a high-activity steady-state op- 
eration, in which the surface temperature is up to 150 
K higher than the nominal inlet temperature. This in- 
creased is the result of the run away heat-up and 
only remains at a finite steady-state value because fur- 
ther heating would lower the surface catalytic activity; 
for the near-ambient conditions shown primarily because 
enhanced CO desorption would deplete the CO popu- 
lation at the surface. This limitation by CO depletion 
also rationalizes why a high-activity branch only exists 
for more CO-rich conditions than the intrinsically most 
active "state" towards the right in Fig. [H and why the 
corresponding steady-state surface temperature progres- 
sively increases with increasing Pqq. As a net result, the 
observable TOF profile shown in Fig. [1] is significantly 
changed from the intrinsic one, with the highest abso- 
lute rates about one order of magnitude higher and at a 
partial pressure ratio that is significantly shifted towards 
more CO-rich conditions. 

The pronounced surface temperature increase is a di- 
rect consequence of the adiabatic boundary at the back 
of the single-crystal and demonstrates the relevance of 
the detailed heat transfer through the substrate on the 
in-situ measurable activity. If we move to the oppo- 
site extreme where an ideal heat coupling is able to 
maintain the nominal temperature at the sample back- 
side, we never obtain any significant increase of even 
when going to higher temperatures around 600 K, where 
the intrinsic TOFs in the most active "state" exceed 
lO'^ site~"'^sec~^.— >£ However, at such high conversions 
an extended product boundary layer rapidly builds up 
above the flat-faced single-crystal surface with concomi- 
tant mass transfer limitations then affecting the observ- 
able TOFs. This is illustrated in Fig. H where at T'°' = 
600 K the observable TOFs are never as high as in the 
intrinsic most active "state" , but with a high-activity 
steady-state branch now extending to much less CO-rich 
conditions. The latter is possible because the boundary 
layer limits the transport of both diatomic reactants in 
a similar way. As long as the nominal inlet composition 
is more CO-rich than stoichiometric feed this effectively 
increases the Pcq/Pq^ ratio directly at the surface as 
shown in the lower panel of Fig. [5] and thereby creates 
an impingement that brings the surface closer to its most 
active "state". 

Even in this ideal isothermal case the true catalytic 
function is thus masked, this time by mass transfer lim- 
itations at the flat-faced surface. In fact, as also shown 
in Fig. [2] the observable TOF is under these conditions 
very close to the upper limit set by pure gas-phase mass 
transfer. This limit is completely independent of the ac- 
tual model catalyst employed, and as apparent from Fig. 
[5] for a wide range of gas-phase conditions it is signifi- 
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FIG. 2: Same as Fig. [1] but for an ideal heat coupling at the 
backside of the crystal. With all other parameters identical 
to those in Fig. [1] the significantly increased intrinsic TOFs 
(black dashed fine) at T'"' = 600K now induce an extended 
product boundary layer above the flat-faced model catalyst 
with the concomitant mass transfer limitations affecting the 
observable TOF profile (red solid line). For a wide range of 
the gas-phase conditions shown the observable TOF profile 
coincides in fact with the upper limit set purely by the gas- 
phase mass transport (blue dotted line), i.e. with a limit 
that is completely independent of the actual model catalyst 
employed. The lower panel shows the corresponding partial 
pressure ratio directly at the surface and illustrates that in 
the high-activity branch, the surface effectively sees a more 
CO-rich gas phase compared to the nominal inlet composition. 



cantly lower than the intrinsic catalytic activity of the 
RuO2(110) surface. With the single-crystal in real ex- 
perimental setups neither perfectly heat-coupled nor iso- 
lated, the two effects discussed here separately in Figs. [T] 
and [2] will obviously be intricately intermingled and need 
to be disentangled by dedicated measurements and se- 
tups. This holds even more as the high-activity branches 
in Figs. [T] and [2] are not the only steady-state solution. 
Instead there is in both cases a range of gas-phase con- 
ditions, in which the system can also operate in a low- 
activity mode that never exceeds the critical TOF and 
therefore coincides with the intrinsic activity. While the 
intrinsic surface kinetics of the employed Ip-kMC model 
does not provide multiple steady-states^, they thus enter 
at larger length scales through the coupling to the sur- 
rounding flow field. A corresponding bi-stability clearly 
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suggests that the system could oscihate between the two 
modes, possibly even inhomogeneously in form of reac- 
tion fronts over the single-crystal surface. In case of the 
first discussed heat transfer limitations, an intuitive prop- 
agation mechanism would hereby be via the formation of 
local hot spots, while in the latter discussed mass trans- 
fer case it would be via gas-phase coupling^, with the 
presented approach establishing the intriguing possibility 
to quantify these model conceptions with first-principles 
based simulations. 



IV. CONCLUSIONS 

We have presented a first-principles based multiscale 
modeling approach to heterogeneous catalysis integrat- 
ing accurate Ip-kMC simulations for the surface kinet- 
ics into CFD modeling of the macro-scale fiow struc- 
tures. While its efficient formulation readily allows ad- 
dressing more complex reactor geometries, we have used 
it here to specifically assess heat and mass transfer ef- 
fects in in-situ studies of single-crystal model catalysts. 
Using RuO2(110) as a representative substrate for the 
frequently studied CO oxidation reaction we find corre- 
sponding limitations to significantly mask the intrinsic 



catalytic function at the high conversion rates reached 
at near- ambient gas-phase conditions. Two crucial and 
hitherto largely unappreciated factors specific to the 
model catalyst reactor geometry are in this respect the 
degree of heat dissipation at the back of the thin sample 
and the propensity to build-up a product boundary layer 
above the flat-faced surface with its high density of active 
sites spread over a macroscopic area. A qualified discus- 
sion of a possibly different surface chemistry across the 
"pressure gap" requires accounting for resulting tempera- 
ture and pressure gradients, both in experimental in-situ 
setups and through predictive-quality integrated models 
as the one presented here. Otherwise wrong mechanistic 
conclusions may be derived which hamper our progress 
towards an atomic-scale understanding of the function of 
heterogeneous catalysts at technologically relevant gas- 
phase conditions. 
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